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In several frustrating systems incommensurate behaviors are often observed. For the S = 1 
bilinear-biquadratic model, we show that the main oscillatory behavior, which is proportional to the 
free edge spins, is eliminated in the incommensurate subphase, considering the average of triplet and 
singlet energy spectra under open boundary conditions. In the same way, the 7r-mode oscillation is 
also removed in the commensurate subphase. Moreover, we find that higher order corrections are 
exponentially decaying from an analysis of small size data. 

PACS numbers: 75.10.Jm, 75.40.Mg, 73.43.Nq 

INTRODUCTION 

Commensurate-incommensurate (C-IC) transitions are interesting phenomena in frustrating quantum spin systems. 
In the Haldane gap systems, incommensurabilities have often been regarded as troublesome problems and have rarely 
been discussed in detail until recently. In inelastic neutron scattering experiments, Xu et o/Q have revealed that a 
quasi-one-dimensional oxide, Y2- x Ca :E BaNi05, has an incommensurate double-peaked structure factor. However, an 
analytical interpretation for incommensurabilities has not been clear. 

The spin-1 bilinear-biquadratic model, 

N 

H = ^hi, hi = Si ■ Sj+i + a(8i ■ S 4+ i) 2 , (1) 
i=i 

which plays a role of prototype of Y 2 BaNi05, has the C-IC change point which is corresponding to the Affleck- 
Kennedy-Lieb-Tasaki (AKLT) point a = an = 1/3 HQ. The AKLT point is solvable and has an energy gap above 
the valence-bond-solid (VBS) ground stated [f|. The VBS state is recently getting more attention with reference to 
quantum entanglements |E H @ • 

Analyzing the structure factor among the commensurate and incommensurate subphases, and also the C-IC change 
point, we have deduced two candidates of the real structure factor unifying these two regions and the AKLT point 
|l(j fill]. Following the S0rensen- Affleck prescription^^, we have constructed Green functions from the candidates, 
and compared them with the energy gap numerically obtained under the open boundary conditions 01 . Thus, we 
have found that the Green function consists of two elements which have anomalies in upper- or lower-half plane. 

In this paper we will show that incommensurate oscillatory behaviors in the incommensurate subphase can be can- 
celed, using triplet and singlet energy spectra. In the same way, we can obtain comparable results in the commensurate 
case. Moreover we discuss higher order corrections considering small size systems. 

SHORT REVIEW ON SINGLET-TRIPLET ENERGY GAP 

The gapped Haldane phase has nonvanishing nonlocal string order an d effectively free S — 1/2 spins at the 
ends of open chains 0, 0|. These edg e sp ins bring a low- lying excitation, triplet, which is degenerate with the 
ground state in the thermodynamic limit 15]. While this excitation is clearly a boundary effect, it is linked to the 
bulk behavior because of the existence of nonlocal string order 0. 

According to S0rensen and Affleck (SA)^], the effective Hamiltonian can be connected with the Green function: 

H cS = (-lfS' L • S'rA 2 J S G («. K ) exp(iqN)S(K), (2) 

where S'l and S'r are spin-1/2 operators at the ends of chain. Considering the Clebsch-Gordan coefficient, we obtain 
singlet and triplet expectation values for edge spins as 

(S|S' L • S' R |S> = -3/4, 

(T|S' L • S' R |T) = 1/4, [ ' 



2 



where |T) = |,s T = 1, s z = ±1,0) and |S) = \s T = 0,s z = 0)[l3- Thus the energy difference between the singlet and 
triplet states can be described with the Green function: 



= (-l) N \ 2 J ^G(q,0)exp(iqN), 

where Et and Eg are the triplet and singlet energies, respectively. 
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FIG. 1: Oscillation of the energy gap of edge states. The energy gaps with different chain lengths (iV = 9, • • • , 13) are plotted 
as a function of a. 

Figure n shows the energy gap behaviors of the model with different chain lengths. Increasing the chain length 
on some fixed a in the commensurate subphase, we see that the energy gap oscillates between even and odd chains 
and decreases exponentially fast, while this behavior is not simple in the incommensurate subphase. The even-odd 
oscillation is modulated by the frustration in the incommensurate subphase. 

In our previous study [13J, which has been performed to explain these behaviors, we have found 



A^sT(iV) 



' (-l) N Aeyyi-mN) sin(VdN), 
0, 

(-l) JV 2exp(-miV) smh(VdN). 



(a > 1/3) 
(a = 1/3) 
(a < 1/3) 



(5) 



The parameters, A, m, and d, depend on a — cud and have been determined with the nonlinear least-squares fitting 
method as shown in Ref . 0] . We have found that the C-IC change is characterized by the following Green function: 



G(q,0) 



1 



1 



(q — im) 2 — d (q + ini) 2 — d 



(6) 



AVERAGE OF TRIPLET-SINGLET ENERGIES 



Now we define the average of triplet-singlet energies: 

£ts = (E s +3E T )/4. 



(7) 



From Eqs. J2J and we expect that we can eliminate the principal term caused by the anomalies of the Green 
function using this average. 

Figure [3 shows bare triplet and singlet energies and averages of triplet-singlet energies plotted as a function of 
l/N for (a) a = 0.3258 (commensurate case), and (b) 0.3425 (incommensurate case). We see that the average of 
triplet-singlet energies on a fixed a varies linearly with 1 /N not only in the incommensurate subphase but also in the 
commensurate subphase. 
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FIG. 2: Cancellation of triplet-singlet energies. Triplet and singlet energy spectra, and an average of triplet-singlet energies 
are plotted as a function of 1/N for 1 < N < 15 when (a) a — 0.3258 (commensurate case), and (b) 0.3425 (incommensurate 
case) . 
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FIG. 3: Finite size effect in an average of triplet-singlet energies and a least-squares fitting line (co + ci/N) are plotted when 
(a) a = 0.3258, (b) 0.3425. 



We perform the least-squares fit to the averages of triplet-singlet energies using the following fitting function: 

f(N) = co + ct/N. (8) 

Figure 13 shows the data of triplet-singlet averages and f(N). When (a) a = 0.3258, each coefficients are obtained 
as c = -0.68167581 ± 2 x 10" 9 and c\ = -1.1671 x 10~ 4 ± 2 x 10~ 8 for 8 < N < 15. In the same way, c = 
-0.6483465831 ± 2 x lO" 10 and ci = -1.77295 x 10~ 4 ± 2 x 10~ 9 when (b) a = 0.3425. The average of triplet-singlet 
energies Ets/N seems to behave high linearly in small a region. In fact, we observe that coefficients of 0(1/N 2 ), 
0(1/ N 3 ) are very small. Hence, we can say that a higher power of 1/N does not appear. 

Figure 01 shows coefficients cq and c\ in Eq. © obtained with the least-squares fitting. The second term in Eq. 
ijSJ comes from surface effects. Since the surface effects are caused by the one-dimensionality, they are different from 
the effect of edge spins. Additionally, the surface effects are perfectly zero, namely c\ = 0, at the AKLT point (lj|. 
Moreover, c\ always shows negative value except the AKLT point. We see that the fitting parameters Co and Ci 
behave approximately as a — «d and (a — arj) 2 , respectively. 

So far we have excluded data smaller than N = 8 since these data differ from Eq. (JSJ. Then we proceed to study 
small size corrections from Eq. (JSJl. For a rough estimate, we consider the difference 

Ay(N) = E TS (N)/N - (c' Q + c'JN) (9) 

for small N, although we use c' and c[, which are determined from the data of N = 14,15. Figure (a) shows 
log(| Ay|) for 1 < N < 8 at a = 0.3258 in the commensurate region. Since the logarithm of Ay decreases linearly with 
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FIG. 4: Coefficients (a) Co and (b) ci in Eq. (SJ are plotted near the AKLT point. It appears that Co oc a ~ an and 
ci (X (a — qd) 2 - 



TV, the difference Ay results in Ay ~ exp(— N/£'). We estimate £' as shown in Fig. [11(b), and then we find £' ~ £/2 
(The correlation length £ is obtained in the previous studv[l^|')[l9|. 
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FIG. 5: (a) The logarithm of |Ay| (Eq. ©) are plotted against TV. The dotted line is drawn for a guide to the eye. (b) Two 
different correlation lengths and the ratio are plotted in the commensurate region. 

In the incommensurate region, we see the oscillatory behavior again. We roughly estimate the wave number of this 
oscillation, and then we find q[ c ~ 2(gic + 7r)|20j. We will show a detailed calculation in another paper. 



SUMMARY 



We have shown that the principal oscillatory behaviors, which are proportional to S'l • S'r, among the triplet 
and singlet energy spectra under open boundary conditions cancel out in the commensurate and incommensurate 
subphases. 

We have found that the energy spectra of singlet and triplet states under open boundary conditions consist of the 
bulk, surface, and edge spin energies: 

E c /N = B + S/N+(C\H cS \C), (C = {S,T}) (10) 

where B and S are the bulk and surface energies, respectively. 

Considering small TV, we have found an exponentially decaying correction term, the correlation length and the 
incommensurate wave number of which differ by factor 2 from those obtained by the energy gap of edge states. One 
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possibility is that the correction will be 0((S'l • S'r) 2 ). Therefore, we will need to improve the SA theory so as to 
contain such higher terms. 

We observe similar results for the S = 1 next-nearest-neighbor model. Thus the cancellation of triplet-singlet 
energies can be found in general spin gap systems. 

The numerical calculation in this work is based on the program packages TITPACK version 2, developed by 
Professor H. Nishimori. 

This research is partially supported by a Grant-in- Aid for Scientific Research (C), 18540376 (2006), from the 
Ministry of Education, Science, Sports and Culture of Japan. 
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